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^ ■ ABSTRACT 

The question concerning the possibility of a first order surface transition in a 
semi-infinite Blume-Capel model is addressed by means of low temperature ex- 
pansions. It is found that such a transition can exist, according to mean field and 
cluster variation approximations, and contrarily to renormalization group results 
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The semi-infinite spin-1 Blume-Capel model has been recently intro- 
duced to describe the surface critical behavior of magnetic systems [1] and 
3 He- 4 He mixtures[2]. The phase diagram of this model has been investigated 
by means of some well-known methods, namely mean field approximation[l- 
3] (MFA), real space renormalization group[4,5] (RG) and cluster variation 
method[6] (CVM). All these methods agree to show that the model can ex- 
hibit three different phases: a completely disordered phase (D), a completely 
ordered phase (O) and a partially ordered phase (S), in which order is local- 
ized on the surface of the lattice. Only very recently it has been recognized [3] 
that, at least at the mean field level, for certain values of the model param- 
eters, two different O phases can be found, with different degrees of order at 
the surface. 

Accordingly, different species of phase transitions are found. Lowering 
the temperature the model can either order as a whole, passing from the D 
to the O phase, and the transition is said to be ordinary, or it can undergo 
two different transitions, a surface one from the D to the S phase and then, 
at a lower temperature, an extraordinary one from the S to the O phase. In 
the limiting case between the two above, the transition is said to be special. 

In the MFA phase diagram[l] all these transitions can be either second or 
first order, depending on the model parameters (only the transition between 
the two O phases considered in [3] is always first order) , and this is confirmed 
by the CVM analysis [6]. On the other hand, the RG scheme developed in [4] 
(but not the one in [5]) excludes the possibility of extraordinary or surface first 
order transitions at finite (i.e. nonzero) temperatures. It is then legitimate 
to ask whether such transitions do really occur or they are only features 
introduced by some particular approximation. 

In the present Letter we have addressed this question by means of low 
temperature expansions. By using a technique similar to that outlined in [7] 
in the context of the Pirogov-Sinai theory (although nothing is rigorous for 
semi-infinite systems) we show that a surface first order phase transition can 
be found at low, but nonzero, temperature, and this also implies the existence 
of an extraordinary first order transition. 
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Our method consists in expanding the free energies of the S and D phases 
taking into account only a few excitations over the corresponding ground 
states. The phase boundary (a first order transition line) is then determined 
by comparison of these free energy. 

Let us consider the Blume-Capel model on a semi-infinite lattice with a 
free surface, which is described by the following hamiltonian: 

-pH = J s ^ SiSj -A S J2 Sf + JbJ2 S ^i ~ A bJ2 S l C 1 ) 

(ij) i (kl) k 

where Si = +1, 0, —1, ^ denotes a sum over all nearest neighbors (n.n.) with 

both sites lying on the surface, denotes a sum over the remaining n.n., and 

(kl) 

f3 = (/c_bT) _1 (with ks Boltzmann constant and T absolute temperature). J s 
and Jb (both positive, since we limit ourselves to the ferromagnetic case) are 
reduced surface and bulk exchange interactions, while A s and Ab are reduced 
surface and bulk anisotropy, respectively. 

In the following we will treat the case J b = J s = J, D = A b / A s > 3/2 on 
a simple cubic lattice with a (100) free surface. In such a case (see e.g. [3]) 
the ground state is O for A b < 3J, S for 3J < A b < 2D J and D for A b > 
2D J and, according to mean field theory, in the (A&/6J, T) phase diagram, 
first order extraordinary and surface transition lines start from (1/2,0) and 
(D/3,0) respectively. Let us now concentrate on the surface transition and 
develop low temperature expansions for the surface contributions and 
fs°^ to the free energies of the S and D phases. 

To begin with, let us consider the S phase, which is characterized by 
Si = +1 (or, equivalently, —1) for all surface sites, and Sk = for all bulk 
sites. On a lattice with iV surface sites, its reduced ground state energy is 
Eq S ^ = N(—2J + A s ). To determine a low temperature expansion for 
up to the 5th order in x = e~ J we have enumerated all the excitations with 
respect to the ground state with energy 5E < 6J. An excitation is regarded 
as a set of spins which are in a state different from the ground state, and 
is represented by a graph, obtained drawing a set of circles containing the 



3 



values of the changed spins and joining with a line each pair of neighboring 
spins. We have then surface, bulk and mixed excitations, depending on which 
spins are changed. To each excitation g (see Tab. I for some examples) we 
have assigned a multiplicity M(g), which is the number of embeddings of 
the corresponding graph in the lattice under consideration and a Boltzmann 
weight e - 5E(9 \ 

We can thus write the partition function in the form 



Z (S) =e -K S) J2M(g)e 



-SE(g) 



(2) 



which, up to the 5th order in x, reads (for a lattice of L layers, each of N 
sites) 

Z^ S \N, L) = 



2w~ N { 1 + N 



2\„3 



(w + y)x z + 2(w 2 + (L — 2)y + y z )x 



+ 



N -5 



-w 2 + 6w s + w 4 + y + 



N -3 



y z + 6y 3 + y 4 + (N - l)wy \x 4 



2 "2 
+ 2^(N- 8)w 3 + 9w 4 + 4w 5 + w 6 

+ (NL -2N + 3L- 8)y 2 + (N - 6)y 3 + 9y 4 + 4y 5 + y 6 
+ (N(L - 2) + l)wy + (N — 2)w 2 y + (N - 2)wy 2 ^)x b J, 

(3) 

where w = e As_2J , y = w~ D x 2e and e = D - 3/2. 

Given a partition function Z{N, L), the surface free energy density of the 
corresponding semi-infinite system is given by [8] 



lim 



In Z (N, L) — NLfb 



where is the bulk free energy density and is given by 

,, = _ lim l ^ N ' L \ 

J N,L^oo NL 

4 



(4) 



(5) 



We then have, for the surface free energy density of the S phase: 



/j s > = A s - 2 J - { (w + y)x 2 + 2(w 2 - 2y + y 2 ) 



x 



+ 
+ 2 



w 2 [ w 2 + 6w - ^ j + (1 - w)y - ^y 2 + 6y 3 + y 4 



x 



w 3 (w 3 + 4w 2 + 9w - 8) 



+ (1 - 2w)wy - 2(w + 4)y 2 - 6y 3 + 9y 4 + 4y 5 + y 6 



x 



In the same way we find 

Z( D \N,L) = 

1 + 2n{w~ 1 



2w~ + (L — l)y 



x" 



+ |w 2 (w 2 + 6w 1 +N-5) + w 1 y 



x 



+ 



2w~ 



and 



(w~ 4 + 4w~ 3 + 9w~ 2 + 2(N - S)^- 1 + l) 
+ 2(NL - N - l)w~ l y + 4w~ 2 y + (3L - 4)y 2 a; 5 } 



f(D ) = _ 2 [w~ x x 2 + (2w~ 2 - y)x 3 

^w~ 2 {w~ 2 + 6uj _1 — 5) + w~ x y 



+ 



+ 



x 



(6) 



(7) 



(8) 



2«r 2 («r 4 + 4w~ 3 + 9w~ 2 - WW' 1 + 1) 

- 2(1 - 2w~ 1 )w~ 1 y - V]x 5 }. 

The transition line is then readily obtained by numerically solving the 
equation = fs°^ for A&/6J at each temperature, and is reported in 
Fig. 1 together with MFA and CVM results, for a choice of values of the 
model parameters. The transition is first order, as can be checked by looking 
at derivatives of the free energies, and the results of MFA and CVM are 
qualitatively confirmed (the reentrance predicted by CVM is not seen here). 
The above results should be particularly reliable for e < 1/12 (in order to 
satisfy SE(g) < 6 J for all the excitations considered) and for /c B T/6J < 0.06, 
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so that the 5th order contributions are no more than a few percent of the sum 
of the lower order terms (excluding the Oth order ones, of course, since they 
are defined up to an additive constant), and these conditions are satisfied in 
Fig. 1. 

Implicitly, since the bulk phase diagram is simply the phase diagram of 
the infinite system, which is known to undergo also a first order transition, 
our results show that the Blume-Capel model on a semi-infinite lattice is 
also capable of exhibiting an extraordinary first order transition. 

Summarizing, we have shown by low temperature expansions that the sur- 
face transition of the semi-infinite Blume-Capel model can be also first order, 
according to previous results obtained by MFA and CVM, and in contrast 
with a RG analysis. We believe that the present method is more accurate 
than the RG scheme proposed in [4] , since this kind of low temperature ex- 
pansion, even if not rigorous in the presence of a free surface, is particularly 
well-suited for first order transitions. Furthermore, RG schemes for models 
with several different phases and rich phase diagrams heavily depend on the 
choice of the mapping truncation, as can be seen by comparing the results 
in [4], which exclude the possibility of a surface or extraordinary first order 
transition, with those in [5] , where, by means of a different RG scheme, such 
transitions are found. Indeed, only very recently [9] Berker and Netz have 
succeeded in devising a RG scheme capable of describing the antiquadrupolar 
and ferrimagnetic phases of the Blume-Emery-Griffiths model. 
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Figure Caption 

Fig. 1: Surface first order transition lines for e = 0.03 as given by the 
present method (solid line), MFA (dashed line) and CVM (dotted 
line) . 
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Table I. Examples of elementary excitations of the S phase 
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